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PROCEDURE FOR NOISE PREDICTION AND OPTIMIZATION OF 


ADVANCED TECHNOLOGY PROPELLERS 

Wen-Huei Jou and Samuel Bernstein 
Flow Research Company 


SUMMARY 

The objectives of this work are to formulate a good method for 

predicting the noise field generated by an advanced technology propeller 

operating at supersonic tip speed and to formulate a noise minimization 

scheme for the design of the propeller blades. The work of Hawkings and 

Lowson on the noise prediction based on the Ffowcs Williams-Hawkings 

equation has been extended to include a foreward flight velocity. The 

main idea for the extension is that the acoustic wave front propagates 

in a uniform flow as a convected and expanding sphere. Based on this, 

the fundamental solution of the convected wave equation can be easily 

obtained. Acoustic pressure at the observer’s position is Fourier 

analyzed to give the frequency spectrum. The Fourier coefficients are 

expressed as integrals over the sources on the blade surface. The 

results show that cones of silence exist fore and aft the rotor plane 

with the semiapex angle ~ .-11 , w 

= sm — , where M, . is the tip Mach 

s M tip tip 

number. This semiapex angle is the same for the foreward and backward 
cones and is independent of the flight Mach number. This result is 
confirmed by the computation of the ray path associated with the emitted 
Mach waves. The strength of the signal is further modified by the 

j 2 2 ^ | —l 

Doppler amplification factor j(l-M sin SI) 2 - M cos and is stronger 

on the downstream side. These explicit directivity features are believed 
to give better physical insight into the problem of noise generation. 

An optimization scheme for the blade design is formulated. The 
practical application of the scheme depends on the ease of the prediction 
of the load distribution from flow field computation. 


INTRODUCTION 

The motivation behind the present interest in advanced propeller 
aircraft is provided by a potential fuel saving on the order of 15 
percent over a turbofan-powered aircraft (Ref. 1). Over the past several 
years, NASA has sponsored several theoretical and experimental studies 
to explore the feasibility of implementing a new propeller design for a 



flight Mach number of 0.8 at 30,000 feet. Several configurations of 8- 
blade (and more recently 10-blade) propellers have been studied, and 
while the required efficiency of these propellers appears achievable, 
the potential noise level associated with these designs is extremely 
high. Estimates of the expected noise level by Hamilton Standard and by 
Boeing (Refs. 2, 3) are on the order of 135 to 150 dB for a tip Mach 
number of 1.2. 

The objective of this present work is to formulate a noise compu- 
tation and optimization procedure to enable the investigation of pro- 
peller design changes to minimize the radiated noise. The analysis was 
directed towards the noise associated with the transonic and supersonic 
flow over the propeller, and it will be interfaced with the flow field 
computations currently underway at NASA Ames Research Center. 

This effort was supported by Professor W. D. Hayes of Princeton 
University who served as a consultant on this project. Professor Hayes 
was primarily responsible for the discussion on wave kinematics. 


Propeller Noise Classification 

The noise generated by a propeller with supersonic tip speed can be 
classified into three categories: 

(1) Multiple pure tone (MPT) : This noise appears with a discrete 

spectrum at the shaft rotating frequency and its multiples, and it is a 
unique feature of a propeller with supersonic tip speed (Ref. 1, 2, 3, 
4). It has been identified as being caused by the merging of the 
rotating shock waves as they propagate away from the blades. Due to 
imperfection of the blade pitch angle and spacing, there will be a 
slight difference in the strength of the shocks as they emerge from the 
blades. When the shocks propagate into the far field, the stronger 
shocks catch up with the weaker ones and they merge into a shock of 
higher strength. Eventually, the shocks merge into a single shock 
rotating at the shaft frequency. MPT is an important component in the 
turbine fan noise. However, with only 8- or 10-blade configurations 
considered for the transonic propeller, this component of noise occurs 
in the very far field where the noise level is sufficiently weak and 
does not cause any concern. Also, it is a component of the noise which 
is hard to eliminate. We shall not pursue MPT in this report. 

(2) Wide band noise: This is noise with a continuous spectrum, 

and it is attributable to turbulence traversing the propeller. For 
supersonic propellers, the wide band noise is much lower than the blade 
passing tone discussed in the following paragraph. 

(3) Blade passing tone (BPT) : This is a discrete tone at the 

blade passing frequency and its multiples. It is the rotating steady- 
pressure distribution sweeping past the observer periodically at the 
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blade passing frequency. With supersonic tip speed, a Mach wave system 
is generated around the blade. The passing of this rotating Mach wave 
system through the observer is the main cause of the pronounced noise 
encountered with supersonic propellers. In the remainder of this report, 
we shall focus the discussion on this component of the noise. 


THEORY OF BLADE PASSING TONE 

The work of Gutin on the noise generated by an open rotor (Ref. 4) 
is the only theoretical work published before the late 1960 f s. The main 
reason is that propellers were then operating at the subsonic tip speed 
range; Gutin f s theory gives a very good description of the noise field 
in subsonic fields. 

With the increasing interest in propellers operating at transonic 
and supersonic tip speeds, Gutin r s theory has been proved to be inadequate. 
As the tip speed approaches and exceeds sonic speed, the characteristic 
wavelength of the noise is no longer much larger than the blade dimensions, 
as implied by Gutin 1 s theory. The concept of representing a blade by a 
rotating concentrated force for calculation of the noise field breaks 
down. The noise field now depends on the detailed distribution of the 
source strength on the blade. 

The general theory of sound generated by a surface in arbitrary 
motion has been investigated by Ffowcs Williams and Hawkings (Ref 5). 

They start from Lighthill's equation for aeroacoustics , with additional 
source terms resulting from the assumption that the solid surface is a 
surface of discontinuity within which the fluid is stationary. The 
general solution can be constructed in terms of the fundamental solution 
of the forced wave equation. The general solution is cast in an integral 
form over both the space and time domain involving a Dirac delta function 
in the integrand. Now, the important concept is that the space-time 
domain can be viewed as a general, four-dimensional space over which the 
Dirac delta function projects the integral to a three-dimensional subspace. 
This three-dimensional subspace can be chosen at will for convenience, 
thus resulting in various forms of representation of the solution. One 
of the most important forms is the projection to the physical space 
domain. 

The procedure results in the traditional retarded potential solution. 
The integrand for integration in the space domain is then shown to be 
singular when the relative velocity of the surface to the observer is at 
sonic speed. One may also project the four-dimensional integral to a 
surface of sphere collapsing at sonic speed toward the observer. This 
form of solution avoids the above mentioned singularity, but creates a 
new singularity when the solid surface is moving toward the field point 
at sonic speed with its normal in the direction of radiation. 

Two major works have been published applying the general theory of 
Ffowcs Williams and Hawkings to the supersonic rotor noise. Hawkings 
and Lowson (Ref. 6) base their analysis on the retarded potential solution. 
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They solve the problem in the frequency space and construct a time-wave 
form from its Fourier components. The singularity in the integrand is 
shown to be integrable. The Fourier coefficients are obtained explicitly 
by using a far field approximation. This leads to an integral over the 
blade surface which is evaluated numerically. 

The authors then go on to apply Whitham’s theory of weak shock 
propagation to calculate the nonlinear distortion of a time-wave form as 
it propagates into the far field. They show that the nonlinear effects 
substantially alter the wave form (hence the spectrum), but have weak 
effects on the overall sound pressure level. Finally, Hawkings and 
Lowson conclude that the linear theory can be used as a prediction tool 
for the gross features of the acoustic field. Since we shall attempt to 
generalize Hawkings and Lowson 1 s theory to include the forward flight, 
we will expand on some features of their theory in the course of our 
analysis . 

Farassat and Brown (Ref. 7), on the other hand, take up the approach 
using the collapsing sphere solution of Ffowcs Williams and Hawkings’s 
general theory. To obtain the time-wave form, a collapsing sphere, 
centered at the observer, is generated from infinity for each time step. 
As the sphere collapses, the intersection of the sphere and the rotating 
blade surface is computed. The intersection line traces a surface, and 
the noise for that time step is evaluated by integration over that 
surface. No far field approximation is made; the procedure is purely 
numerical. The noise field can be mapped out by calculating noise at 
different observers’ positions. 

For the purpose of optimization, we are more interested in a scheme 
which contains the following features: 

(1) Fast computation of the noise field even if some approximation 

must be made. 

(2) Some explicit directivity distribution as a guidance for 

choosing observer’s position for noise minimization. 

(3) Correct prediction of the sensitivity of the noise field to 

propeller design parameters. 

With this in mind, we have elected to use Hawkings and Lowson f s approach 
in outlining a noise minimization procedure. 

In the following section, we extend the work by Hawkings and Lowson 
to include a forward flight. The resulting analysis will be applied to 
compute the noise level at the aircraft cabin, which is one of the major 
concerns for the advanced propeller aircraft. 
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PROPELLER NOISE WITH FORWARD FLIGHT 
Retarded Potential Solution of Convected Wave Equation 

The basic idea in the following analysis is that the wave front 
propagates in the uniform flow as a convected sphere (Ref. 8). Based on 
this, we are able to construct fairly easily the fundamental solution of 
the convected wave equation and cast the general solution in an integral 
form. 


Consider an acoustic disturbance generated at A and propagating 
in a medium with speed of sound a and convecting Mach number M . 

The signal propagates away from A as a convected sphere, as shown in 
Figure 1. 



Figure 1 . Wave Front Geometry. 


The center of the convected spherical wave front is at B . For an 
observer at C , the received strength of the signal is proportional to 
f-1 rather than r~* . By simple geometry, the relation between r and 
r can be obtained as 


5 


K 


( 1 ) 



(K 2 + B 2 )* 5 - 


where 


with 


The unit vector pointing 
time to the observer at 


K = (M • e r )B Z 

3 = (1 - M 2 ) -1 
2 

< 1 . 

from A toward C is 
C is 


( 2 ) 

(3) 

(4) 

e . The retarded 
r 


T 

O 



(5) 


The starting point of the subsequent analysis is that the solution 
to a convected Ffowcs Williams-Hawkings-equation, 


1 D 2 „2 3 n 

2 p p 3x. 3x.... ^ij 


a Dt 


i 3 


is the following: 


4Tr(p - P o ) = 3xjx .... / Q ij 


A(t - t q ) ^ 

• • • (y , t ) dy dT 


with the relations (1) ~ (5), 


( 6 ) 


(7) 


In the preceding solution, the source space y is described in a 
coordinate system fixed in space. The moving body (emitter) is viewed 
as distributed sources over the entire space and their strength is non- 
zero only when the solid body passes through that field point ($) . 

In other words, the sources are ,f pulsing sources” distributed in the 
entire source space. 


In the case of the acoustic emission due to the thickness of and 
forces on the blade, the source space is chosen best in a coordinate 
fixed in the blade. The source strength described in such a coordinate 
system will be in a steady state. The acoustic signal at a point in the 
observer’s space 5 and at the observer’s time t can then be obtained 
by summing the contributions of the steady sources emitted at the position 
related to the observer’s space and time (£, t) by the retarded time. 
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To do this, a transformation of the source space from the inertial 
frame y to the blade fixed coordinate is required- The transfor- 
mation is given by: 


V2 



Where 


y 2 ” ^3 sin a3T + ^2 COS WT ’ 

y 3 * ^3 COS WT ” ^2 sin WT * 
and T * t — — | (K 2 + B 2 )* 5 - K 


Without loss of generality, we assume the initial phase of the rotation 
is zero- Since r and e depend on the variable y , the above 
transformation is implicit- By some tedius calculations, it can be 
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shown that the Jacobian of the transformation is: 


J = det 


9y i 


1 


(1 + u i 

J 1 


( 8 ) 


where U-£ is the velocity of the emitter. Note that for M=0 , so 
that K=0 , and 3=1 , we recover the result shown by Ffowcs Williams and 
Hawkings (Ref. 5) and by Lowson (Ref. 9); i.e.. 


J - 


1 


The acoustic pressure is then 


_1 

U + 
— * e 
a r 


8 f 

1± j 

dS ». f 

p^a(m + M) • Vh 

3x. I 


dS St / 

J 

A 

3 J 

r 

J 

r 


(9) 


, ( 10 ) 


where denotes the values of the integrand at the retarded time, 

and the integration is evaluated over the blade surface. A thin blade 
is assumed, and we sum the force on the top and bottom surfaces of the 
blade. The net force on the planform of the blade is f j ; the thickness 
of the blade is h ; and the rotation Mach number of the emitting element 
is i . 


Fourier Components of Acoustic Pressure 

For the reasons mentioned in the Section "Theory of Blade Passing 
Tone" of this report, we follow Hawkings and Lowson ? s analysis (Ref. 6) 
and expand the solution in Fourier series. To convert the Fourier 
integral from being over the observer's time domain to the emitter’s 
time domain, we have 


so that 


T + At (y(x)) = t , 


dtt ^liM iT , it 

dT ay ± 


(ID 
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I 


But, 


so that 


or 


dy i „ , 3 (At) 3 t 

- t — = U. and -5 = - 5 — , 

dx 1 3y 1 3y i 


dx 


( 1 + 0 i%)- 


dt 


dt = J 1 dT 


( 12 ) 


With this, the Fourier coefficients can be shown as 


C 

n 


j__ 

4tt 


J (intol 


9 


9x . 

3 


!_.) dS 


(13) 


where to is the rotation frequency, n is the harmonic number, S 
is the blade surface. 


I 


to 

TT 



p a(m 4- M) 


Vh 

— exp 



(T + 


*) 


dT 


(14) 


and 


2t r 

to f. / ^ \ 

x j • ? / 7 1 exp \ int0 (T + f) j dT 

o 


(15) 


9 


I 



3 

For a propeller with supersonic tip speed, the wave length ^ is 
comparable to the length of the blade R , since ^ - 0(1) . The 
source is noncompact, and it cannot be considered as a concentrated 
source even in the far field. However, we consider the far field 
approximation as the limit R/r « 1 and << 1 . The approximate 

expression for r up to the first order is needed. 

The acoustic field is axisymmetr^c and we shall consider an observer 
on x-y plane with position vector r Q . Referring to Figure 3, we 
have 
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So the exponential factor in equation (14) can be approximated by 


with 


and 


Let 


with 


E = exp ( iconT + 


^ iconT + icon ~ ^ 


exp 


moo 


(1 


— 2 ~~ F ( 1 - M 2 sin^ 
- M Z ) L 


ft) - M cos ft 


• exp / iconT 


( . inu)R sin ft cos 0 ) 

ia)nT a n w 2 • 2 ( * 

(1 - M sxn ft) ) 


0 = COT + <j> 


sm ft = 

r 

o 


coR 


= m be the rotation Mach number of the emitting element. 


Y 


m sin ft 

(1 - sin^ ft) 2 


and the Doppler factor 



We have* 


2p a (m + M) • Vh (1 - M ) 
o 

r o *0 


incor 


exp 


2-J Kj - in* - [ J n (nY) 


a(l - M ) 


where J is the Bessel function. For the noise generated by the 
n 

forces, we have 


31, 

_ co 

3x. " 7T 
3 

^ 2t t 


. 2 tt 

(0 


f . 


j . 

3 


~ exp ^ inco (t + —■) j 


dT 


(0 

TT 


03 


O 


/ . , ^ rA . 3r / icon 1 

exp ( (T + ; >) 


dx . 


(17) 


For the far field « 1 , we can neglect the last term in the last 

no3r 

2tt 


expression: 


31 . .2 

inco 


3x . aTT 
3 


03 


( A V f . 

iwn (T + f) j j 1 3^ 


dx 


e 


It can be shown that 

1 3r - 1 ( Ig 

r 9x. K^r ) ^2 . 2 

j Do ( (1 - M sm ft) 

The force on the blade is 


M sin ft 


(1 - sin^ ft)^ 


TT J " Mi 


f = (-T , -D sin 0, D cos 0) , 


t r 

* ~ | exp £in - inY sin £ d£ = J^CnY) 

o 
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where T and D are respectively the thrust and drag components, so 
that 

f • e = - T cos ft - D sin ft sin 0 , 

r 

f » Vr _ 1 

■r O O U 

K^r^ (1 - M sin ft) 2 

|t J^- cos fl + M( 1 - M 2 Sin 2 J - D sin ft sin 0 (1 - M 2 ) | . 


We have 


31, 


3x . 
J 


ina) 

aTT 


r 2tt 

OJ 


exp ( ino) (t + — ) ) 

' 3 ' K_r (1 


V 


o o 

M sin ft) 2 


!■[- 


2 2 1 
cos ft + M( 1 - M sin ft) 2 J 


- D (1 - M ) sin ft sin (tOT + <)>)> dT 


The thrust component is 


31 . 

,°Vt 


_ 2iruu [- cos ft + M (1 - M 2 sin 2 ft) 2 ] 

3 K D r Q (1 - M 2 sin 2 ft)* 5 


inur 


inTT 


• exp { ~ Kp - in$ - ± y~ > J (nY) 

1 a(l - IT) ‘ n 


„ . .xncor 

a ^ n K exp{ — — 2 — in(J> - (nY) 

o' 


|a(l - M 2 ) 

cos ft 

(1 - M 2 Sin 2 ft)* 5 


inTTl . 

2 r 


- M > T 


For the drag component we have 


(18) 
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/M 2inu 

\ / D “ 


- 2 into D(1 - M ) sin fl 


,2 .2 


exp 


K^r o (1 - M sin ft) 
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( lna)K D r o 
| a(l - M 2 ) 


- in<£ - 


inTT 


cos (nC - nY sin C) cos C dC 


But 


/' 


cos C cos (nC - nY sin 5) dC 


/ TT 

| cos Jj 


(n + 1)C - nY sin C 


+ cos (n - DC - nY sin C dC 


i" «]| 


■f j j „ + i (nY > + j „-i ( " Y > 


So , 


- T J n (n ? ) . 


/ 91 j \ _ 2inoJ D(1 - M 2 ) 
\ 9x j / D = aK D r o m 


exp 


ina)K p r o 

a(l - M 2 ) 


- incj) - 


xnTT 


• J n (nY) . (19) 

Substituting equations (17), (18) and (19) into equation (13), we have 


„ mo) 

C = 0 exp 

n 27taK D r o 


xntor 


oS 


inTT 


2 2 
a(l - JT) 


•Jj |p o a 2 (m + M) • Vh (1 - M 2 ) + T 


cos Q 


9 9 U 

(1 - M z sin Z «)* 


- M 


2 \ 

- = - > J (nY) exp (-in4>) (cos 0) ^ R dR d(J> 

min 


( 20 ) 
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where m, T, D and h are functions of the coordinates (4>,R) on the 
blade surface, and 0 is the twist angle. 

We have reduced the Fourier coefficient of the noise at a point 
x to a double integral on the blade planform surface. 

One interesting feature of equation (20) can be observed. For fixed 
R , the <f) integral is the Fourier integral of the thickness h and the 
forces T and D . With proper arrangement of the planform sweep, the 
R integration may give a favorable phase cancellation effect for those 
particular harmonics. 


Zones of Relative Silence 

Some of the directivity of the acoustic disturbance can be extracted 
from equation (20). In particular, the Bessel function J n (nY) is expo- 
nentially small for Y < 1 . This is especially so for large n . So, 

Y <. 1 defines a zone of relative silence (see Figure 4) in which only 
the subsonic Gutin noise exists: 

2 

y 2 _ m sin SI 

Y 2 2 — 1 • 

1 - M sin « 


or 


n < 


sm 


-1 


(m + M ) 




sin 


-1 




Figure 4. Zone of Silence. 


( 21 ) 
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This result is in agreement with the following kinematic consideration 


At each point on the supersonic portion of the blade, one can find 
the ray path of the acoustic waves. 


m 



Figure 5. Ray Path. 




I 



We obtain thereby 

e i = 0 2 = sin -1 — - 1 ^ (25) 

(m 2 + M 2 ) 

From equation (20) , it can also be seen that the strength of the signal 
decays as (Kj)Tq)-1 , where Kj) is the Doppler factor which is dependent 
on the direction ft as 

^ = (1 - M 2 sin 2 ft)* 5 - M cos ft (26) 

So* the strength of the signal is stronger on the downstream side as one 
would expect intuitively. 
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Figure 6 shows the qualitative directivity distribution from the 
forgoing discussion. A similar picture is obtained using kinematic 
consideration by Lowson and Jupe (Ref. 10). 
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Asymptotic Evaluation of Noise Field for Large Wave Number 

Since we are expecting a spectrum rich in high harmonics, an approximate 
solution may be taken for the high frequency components as h -*■ 00 . Following 
again Hawkings and Lowson (Ref. 6), we use the stationary phase method. 

The phase function under consideration is 


f = U)(T + f> = e + ^ , 

a a 


(27) 


where r is a function of (R,0), and 


A. ~ 

r ~ 


(1 - M Z ) 


*D 


(1 - M 2 ) sin 0 R 
(1 - M 2 sin 2 ft) 1 * r o 


cos 0 


) 
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The stationary point (actually, a line) is: 


n _ _3tt toR = (1 - M 2 Sin 2 ft)** 

2 5 a sin ft 

Cjl)R 

Since — = m , one can see from equation (28) that the main 
contribution to the integral comes from a constant radius line where a 
ray is pointing directly toward the observer, as shown in Figure 7. 


Rays Not Pointing 



i 

I 

Rotor Plane 

Figure 7. Radius of Stationary Phase and Ray. 


For the thickness component, we obtain 



(1 - M sin Q)^ (1 - M ) 


sin f2 
(m + M) • Vh 


Vo 


exp 


i3n7T 


R 

cos 6 Jr=R 


exp (-in<{>) d<f> 


+ 


ina3r o K P 
a(l - M 2 ) 


with 


R 

o 


a (1 - M 2 sin 2 ft)** 
to sin ft 


(29) 
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For the force components, we have 


= - ~^2 JJ R dR d< f> J d< p f □ ^ Hr [^f) ex p < e + f> 


_ i__ (1 - sin 4 * Pi)'* 

2tt sin PI exp 


. 0 int)r K_ 
i3n7T o 1) 


2 ' 2 
a(a - M ) 


f i ! 

1, 

- COS ^ 1 M 

d* 

i-i 

O 

( 

2 2 i- + M 

(1 - M sin ft) 2 


D sin g (1 - M j R ( -in<t ,, 

~ 2 n ,Ss ^Tq{ e d(j) 

(1 - M sxn Pu) ) 

' R— R 


(30) 


Combining equations (29) and (30) results in the following: 


„ _ i (1 - M 2 Sin 2 Sip (1 - M 2 ) 

^ ~ ^ ~ exp 
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"n 2 t t sin ft 
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inTT o TD 


2 ’ 2 
a (1 - M Z ) 


2 -> T 

P Q a (m + M) • Vh + — 




COS P 


.2, l ,, .,2 . 2 " M 


(1 - M ) \(1 - M sin .ft)' 


D sin Pi 


,2 2 


(1 - M sin fi) 2 J 


R ( -inc() , , 
coTe ( e dt f> 
R=R 


(31) 


The asymptotic evaluation of the Fourier integral reduces the surface 
integral to a line integral on the chord, from where the main contribution 
comes. This reduces the computing time needed for the high frequency 
components of the spectrum. 
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OPTIMIZATION PROCEDURE FOR NOISE REDUCTION 


Numerical optimization techniques have attained a high degree of 
refinement. However, very little is known to aid the selection of a 
particular optimization scheme for a given problem. In this study, we 
have evaluated several procedures to minimize the propeller-radiated noise 
which in principle could be implemented with any numerical optimization 
algorithm. The evaluation, however, was guided by our experience of 
utilizing the method of feasible directions and the conjugate gradient 
technique, as developed by Vanderplaats (Ref. 11) for constrained minimi- 
zation problems. In this algorithm, the numerical search procedure is 
performed in the gradient directions to minimize the objective function 
while satisfying a set of constraints. For the advanced propeller problem 
we wish to minimize the noise level at a given point, 

oo 

P(x,t) = y^e~ lnaJt C n (aj,x) , (32) 

n =-co 

by varying some design parameters under some constraints. We assume that 
p(x,t) can be minimized by minimizing the Fourier components C n (u)) , 
and we will focus the discussion on minimization of C n (u)) for a single 
harmonic frequency n . We have identified several procedures applicable for 
advanced transonic propellers as follows: a general noise minimization 

procedure involving all design parameters; an approximate noise minimization 
with propeller sweepback; and noise minimization for high frequency 
noise. In these procedures we assume that the observer is moving with 
the source (i.e., the observer is in the aircraft cabin) at a fixed 
lateral distance shortly aft of the propeller. Such a location is expected 
to have the highest radiated noise level, as illustrated in Figure 6. 

General Noise Minimization Procedure 

The noise function is given in terms of its Fourier components 
C n (o)) and expressed in equation (20) which is repeated here for con- 
venience: 
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C = ^ — exp 

n 27TaK D r o 


moor 


0*0 


mTT 


2 2 

a(l - M Z ) Z 


• I ( i p a^ (m + M) • 

N2 * 

JJ v Thickness 


Vh(l - M ) 


+ T 


cos ft 


(1 - M sin ft) 




- HQ | j ( n y) exp (-in(f>) (cos 0) ^ R dR d<J> 

m i n 


✓ 's 


Thrust 


Drag 


( 20 ) 


As identified, the first term in the surface integral of the equation 
includes the thickness distribution h , which is a function of its 
radial and angular location R and $ (see Figure 3). Likewise, the 
next two terms are the thrust and drag contributions. The incremental 
thrust and drag are also functions of each radial and angular location; i.e.. 


h = h(R,<j>) 
T « T(R,<f>) 
D = D(R,(f)) 

and finally the twist or pitch angle 0 
radial station R : 

0 = 0(R) 


, (33a) 

, (33b) 

, (33c) 

is a function of the propeller 

(34) 

angular stations, the 


If we describe the propeller with N radial and M 
optimization problem can then be stated as follows: 


minimize C (u)) by varying the N(3M + 1) 
n 

design variables h(R,()>) , T(R,4>) , D(R,(f)) , and 
0(R) subject to the following N(M + 2) + 2 constraints: 
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(35a) 


h(R Q ,X) h c (R o »X) (where X is the blade sweep) , 
J T ( R q » <}> ) d<t> < J T c (R o ,4>) d4> , 

ff T (R,cf>) R dR dcf) >_ T_ 


and 


ff R 2 d(R, 


— O 


0) dR d</> £ Q o 


(35b) 

(35c) 

(35d) 


The first two N constraints express the structural requirement of 
minimum thickness and, having the loading below either a structural or 
aerodynamic maximum value at each radial location, the final two integral 
constraints express the integrated thrust and torque values to achieve 
the desired propeller performance. 


The practicality of solving this problem will largely depend on the 
time required to compute the flow field around the propeller and to compute 
the gradient of the flow variables with respect to the N(3M + 1) design 
variables. The simpler problem obviously is to consider only the variation 
of the propeller planform to noise as expressed by the thickness distribution. 


Approximate Noise Minimization Procedure 

In this approach we consider only the contribution of the thickness 
to the noise (i.e., only the first term in equation (20)). A first 
order approximation of the planform change by sweeping the propeller 
backward could be treated by computing the noise without regard to the 
coupled changes in the thrust and drag. This simplification essentially 
eliminates the flow field computation of thrust and drag and may provide 
a more feasible optimization procedure. In this procedure, the objective 
function to be minimized includes only the first term in the integral 
and only MN design variables are considered: h(R ,cf>) and Q(R) . 

The constraints are also reduced to N constraints on minimum thickness 
distribution at each radial location. 

Another approach for an approximate optimization scheme may be 
considered by using the asymptotic solution to the high-frequency noise, 
as illustrated by equation (26). 


Noise Minimization for High-Frequency Noise 

In the asymptotic solution for the high-frequency noise (see equation 
(31)) most of the acoustic signature is assumed to be from a ray directed 
at the observer. The surface integral is then reduced to a line integral 
at the location for which the ray is directed towards the observer. Either 
the general noise minimization scheme or the approximate solution may 
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be applied with the asymptotic noise solution with considerable simpli- 
fication. It should be recognized, however, that the solution would be 
a local solution relative to one observer, and numerous locations should 
be considered to redesign the overall propeller for optimum noise. 


CONCLUSIONS 

The prediction scheme outlined in this report is an extension of 
Hawkings and Lowson f s (Ref. 6) work and is intended to facilitate the 
basic tool for the noise optimization in the design of advanced propellers. 

The method gives some explicit directivity distribution. It gives 
good physical insight into the noise generation, particularly in defining 
the existence of the cone of silence, the Doppler amplification factor 
K d and the critical radius for high-frequency components. In terms of 
numerical computation, the procedure is straightforward and fast, a 
factor of considerable importance for numerical optimization. 

The method, however, is limited in some aspects. Except for the 
work by Hawkings and Lowson (Ref. 6), all the existing prediction methods, 
including this work, are based on linear theory. The nonlinear distortion 
of the wave form as it propagates from the emitter is not included. As 
pointed out by Hawkings and Lowson, the nonlinear effect is important in 
determining the spectrum, while linear theory may be sufficient in 
predicting the overall pressure level. The simplicity of the technique, 
however, is quite attractive for design purposes. A second limitation 
is due to the employment of the far field approximation. The observer 
at the cabin wall is roughly at 2.6 times the radius of the propeller. 

At such a close proximity the higher order term in the far field expansion 
may contribute significantly to the noise. Inclusion of the higher 
order terms can be carried out in the future. In many cases, the first 
term of an asymptotic expansion gives fairly good results for moderate 
parameter values. For the optimization scheme, the important feature is 
the ability to predict correctly the sensitivity of the function to the 
various design parameters. It is recommended that the initial optimi- 
zation scheme be applied to the simplest noise prediction scheme. 

A presentation of this work has been included in the Appendix. 


Flow Research Company 

A Division of Flow Industries, Inc. 
Kent, Washington 98031 
December 1, 1978 
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THE APPENDIX 
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FLOW RESEARCH PRESENTATION NO. 104 


PROCEDURE FOR NOISE PREDICTION AND 
OPTIMIZATION OF ADVANCED 
TECHNOLOGY PROPELLERS 


Wen-Huei Jou 
and 

Sam Bernstein 



OBJECTIVES OF CONTRACT 


Review Existing Theories of Propeller Noise 
Formulate a Prediction Scheme 


Formulate an Optimization Scheme for Noise Reduction 



OUTLINE OF PRESENTATION 


Introduction: Past Analyses of Propeller Noise 

Extension of Hawkings and Lowson's Work to Include Forward 
Flight 

Optimization Procedure for Noise Reduction 


Conclusions 


CLASSIFICATION OF PROPELLER NOISE 


• Multiple Pure Tone 

• Wide Band Noise 

• Blade Passing T one 
Gutin's Theory: 

• Rotating dipole at center of force 

• Theory breaks down for supersonic speed 

XIR = (a/co)/R = a/(coR) > 1 



GENERAL THEORY OF FFOWCS WILLIAMS 

AND HAWKINGS 


Sources 


Volume Displacement Monopole 

Surface Forces Distributed Dipoles 

Stress Field Distributed Quadrupoles 



GENERAL SOLUTION 
(FFOWCS WILLIAMS AND HAWKINGS) 

4it(p-p 0 ) = d n /3xj3xj... J Qjj (y,r)[(5(T-T 0 )/r]dy dr 

r = |x-y| 
t 0 = t-r/a 

• General Space-Time Integral 

• Dirac Delta Function Projects the Integral to Three-Dimensional 
Subspace 




• Subspace Can Be Chosen at Will 



GENERAL SOLUTION (CONT.) 


Choose Physical Space as Subspace - Retarded Potential 

Choose a Sphere Centered at Observer and Collapsing at Sonic 
Speed (Sphere of Coincidental Arrival) 

- I Surface Solution or Collapsing Sphere Solution 



HAWKINGS AND LOWSON 




• Method 

• Retarded potential 

• Blade fixed coordinates for source space 

• Fourier analysis 

• Far field approximation 

• Results 

• Qualitative directivity 

• Nonlinear distortion of time-wave form 

• For design purposes, linear theory is sufficient as a 
prediction tool 



FARASSAT 


Method 

• Collapsing sphere 

• Time domain 

• No far field approximation 

Results 

• Complete numerical computation for every field point 

• No qualitative directivity 



PROPELLER NOISE WITH FORWARD FLIGHT 


• Uniformly Convecting Velocity 

|l /a^D^/Dt^)- p= 0 n /0Xj0Xj... Qjj 

• Wave Front is a Convected Sphere 



K = (M • e r ) p 2 
p = ( 1 - M 2 )' 1 


-<] 



GENERAL SOLUTION TO CONVECTED 

WAVE EQUATION 

4n(p-p 0 ) = 3 n /3xjdxj... jQjj (y,r)[d(T-T 0 )/f]dydT 



TRANSFORMATION OF SOURCE SPACE TO 
BLADE FIXED COORDINATES 



RETARDED POTENTIAL SOLUTION ON 
BLADE FIXED COORDINATE 

4n (p - p 0 ) = -a/3xj|[fj/fj]dS 

-did t j [e 0 a(fri + M) • Vh It ■ J]dS 


S: Blade Surface 



FOURIER COMPONENTS 


C n = / (p - p 0 )e' n<;ot dt 

= J Ip- p 0 )e inc ^ T + r ^ a '(dt/dr) dr 
t + At tylr)] = t 
dt = J"^ dr 


-r 

H 



I 


ro 


FOURIER COMPONENTS (CONT.) 


• Polar Coordinate 



• Far Field Approximation 


Wx,y,o) 


(R/r) « 1 ; (a/ncor) « 1 


• C n - <C n >h + (C n )f 



FOURIER COMPONENTS (CONT.) 


<C n )h = A/(K D r 0 ) J j (fn + M) • VhJ n (nY) exp (-inflfR/cos ©jd+dR 


Y = msinQ/(1 -M^sin^Q)^ 



r 


4 =“ 

-r 


QUALITATIVE DIRECTIVITY 


I. Zone of Relative Silence 


J n (nY) 



Q < sin' 1 [l/(m 2 + M 2 ) % ] = sin' 1 [l /M tip 



QUALITATIVE DIRECTIVITY (CONT.) 


Wave Kinematics (Hayes) 


m 



6 -j =02 = sin'll /(m^ + M^) 1 / 2 ] 




QUALITATIVE DIRECTIVITY (CONT.) 


Wave Kinematics (Cont.) 



QUALITATIVE DIRECTIVITY (COIMT.) 

II. Doppler Amplification 

1/Kp = j (1 - sin^ Q) 1 / 2 - M cos Q } 



*r 

*-3 



HIGH FREQUENCY APPROXIMATION 


§ 


R 0 = a/co[(1 - M 2 sin 2 Q) Vl /sin Q 


Rays Not Pointing 



Rotor Plane 


SUMMARY OF NOISE COMPUTATION 

PROCEDURE 

00 

P(x,t) = X) e inojt C n (co,x) 
n = -oo n 


C n (co) = K(co, K d , M, n )fj je 0 a 2 (m + M) • Vh(1 - M 2 ) 

Thickness 

+ T [cos <2/(1 - M 2 sin 2 Q)' A - M] - D(1 - M 2 )/m 


Thrust 


Drag 


• J n (nY) exp (-in<|>) (R/cos 0) dR d<|> 


h = h(R,<(>) 
T = KRi) 
D = D(R,4) 
© = 0(R) 


OPTIMIZATION PROCEDURE FOR 
NOISE REDUCTION 


Minimize C n (a>) by Varying h(R,<|>), T(R,<|>), D(R,<|>),and 0(R) 
Subject to the Following Constraints: 

• Thickness distribution is longer than minimum value at 
each radial location (structural) 

• Integrated thrust at each radial location is lower than a 
given maximum allowable load (structural/aerodynamic) 

• Integrated thrust sufficient for aircraft propulsion 

• Integrated torque compatible with available engine power 
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